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DRIVER WAVEFORM MODELING WITH MULTIPLE EFFECTIVE 
CAPACITANCES 

BACKGROUND OF THE INVENTION 
The present invention relates generally to 
electronic circuit design using computer simulation 
techniques. More specifically, but without limitation 
thereto, the present invention relates to modeling the 
output waveform of a circuit element driving a distributed 
resistance-capacitance network. 

Methods for modeling the output waveform of a 
circuit element or cell driving a resistance-capacitance 
network generally use a Thevenin equivalent model, i.e., a 
voltage source in series with a resistance, and a load 
having a single effective capacitance. However, previous 
methods produced unsatisfactory results using a computer 
circuit simulation program such as SPICE because the 
actual waveform delays were found to be more than ten 
percent larger than those predicted by the simulation 
program. A method is therefore needed for modeling the 
output waveform of a cell driving a resistance-capacitance 
network that compares more closely to the actual output 
waveform. 

SUMMARY OF THE INVENTION 
The present invention advantageously addresses 
the needs above as well as other needs by providing a 
method for modeling the output waveform of a cell driving 



a resistance-capacitance network that includes multiple 
effective capacitances. 

In one embodiment, the invention may be 
characterized as a method of calculating Thevenin 
5 parameters that includes the steps of (a) initializing 
estimates of effective capacitances Ceffl and Ceff2, of a 
switching threshold delay to, and of a slope delay del tat; 
(b) solving ramp response equations for tO and deltat as a 
function of Ceffl and Ceff2; (c) comparing the estimates 
10 of to and deltat with solutions for tO and deltat found in 
step (b) ; and (d) replacing the estimates of to and deltat 
with the solutions for tO and deltat if the solutions for 
to and deltat have not converged to the estimates of tO 
and deltat. 

15 In another embodiment, the invention may be 

characterized as a computer program product that includes 
a medium for embodying a computer program for input to a 
computer and a computer program embodied in the medium for 
causing the computer to perform at least one of the 

20 following functions: 

(a) initializing estimates of effective 
capacitances Ceffl and Ceff2, of a switching threshold 
delay tO, and of a slope delay deltat; 

(b) solving ramp response equations for tO and 
25 deltat as a function of Ceffl and Ceff2; 

(c) comparing the estimates of tO and deltat 
with solutions for to and deltat found in step (b) ; 

(d) replacing the estimates of tO and deltat 
with the solutions for to and deltat if the solutions for 

3 0 to and deltat have not converged to the estimates of tO 
and deltat; 
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(e) repeating steps (b) , (c) , and (d) until the 
solutions for to and deltat converge to the estimates of 
to and deltat; 

if] calculating a delayl as a function of 
5 t30(Ceffl) or t70(Ceffl) and a delay2 as a function of 
t50 (Ceff2) from a Foster or a pi model; 

(g) comparing delayl and delay2 to delays 
delayl ' and delay2 ' corresponding to Ceffl and Ceff2 in a 
delay lookup table ; 
10 (h) finding new values for Ceffl and Ceff2 from 

a reverse lookup of delayl and delay2 in the delay lookup 
table if the calculated values for delayl and delay2 have 
not converged to delays delayl ' and delay2 ' ; 

(i) replacing the estimates of Ceffl and Ceff2 
15 in step (b) with the new values for Ceffl and Ceff2; and 
(j) repeating steps (b) through (i) until the 
calculated values for delayl and delay2 converge to delays 
delayl ' and delay2 ' . 

The features and advantages summarized above in 
2 0 addition to other aspects of the present invention will 
become more apparent from the description, presented in 
conjunction with the following drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

2 5 The above and other aspects, features and 

advantages of the present invention will be more apparent 
from the following more specific description thereof, 
presented in conjunction with the following drawings 
wherein: 

3 0 FIG. 1 is a schematic diagram of a Thevenin 

equivalent model of the prior art; 
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FIG. 2 is a schematic diagram of a Foster load 
synthesis model of the prior art for modeling the voltage 
response of the CMOS gate of FIG. 1; 

FIG. 3 is a schematic diagram of a pi load 
synthesis model of the prior art for modeling the voltage 
response of the CMOS gate of FIG. 1; 

FIG. 4 is an ideal voltage plot of the prior art 
illustrating Thevenin voltage parameter definitions for 
the CMOS gate of FIG. 1; 

FIG. 5 is a voltage plot illustrating typical 
output voltage response of the CMOS gate of FIG. 1 as a 
function of time; and 

FIG. 6 is a flowchart of an embodiment of the 
present invention for calculating Thevenin parameters and 
the effective capacitances Ceffl and Ceff2. 

Corresponding reference characters indicate 
corresponding elements throughout the several views of the 
drawings . 

DETAILED DESCRIPTION OF THE DRAWINGS 
The following description is presented to 
disclose the currently known best mode for making and 
using the present invention. The scope of the invention 
is defined by the claims. 

FIG. 1 is a schematic diagram of a Thevenin 
equivalent model 100 of the prior art. Shown are a 
Thevenin voltage source 102, a driver resistance 104, a 
CMOS gate 150 and a load impedance 152. The CMOS gate 150 
includes the Thevenin voltage source 102, which is 
connected in series with the driver resistance 104. The 



Thevenin equivalent model 100 is used to model a typical 
cell in integrated circuit designs. 

FIG. 2 is a schematic diagram of a Foster load 
synthesis model 200 of the prior art for modeling the 
voltage response of the CMOS gate 150 of FIG. 1. Shown 
are a driver resistance 104, a first current branch 2 02, a 
first resistor 204, a first capacitor 206, a second 
current branch 208, a second resistor 210, and a second 
capacitor 212. The Foster load synthesis model 200 is 
connected in series with the driver resistance 104 of the 
CMOS gate 150 and includes multiple current branches 2 02 
and 208 connected in parallel. In the first current 
branch 202, the first resistor 204 is connected in series 
with the first capacitor 206. Likewise in the second 
current branch 208, the second resistor 210 is connected 
in series with the second capacitor 212. 

FIG. 3 is a schematic diagram of a pi load 
synthesis model 300 of the prior art for modeling the 
voltage response of the CMOS gate 15 0 of FIG. 1. Shown 
are a driver resistance 104, a first capacitor 3 02, a 
first resistor 304, and a second capacitor 306. The pi 
load synthesis model 3 00 is connected in series with the 
driver resistance 104 of the CMOS gate 150. The first 
capacitor 3 02 is connected between the driver resistance 
104 and ground. The first resistor 3 04 and the second 
capacitor 3 06 are connected in series across the first 
capacitor 302 . 

FIG. 4 is an ideal voltage plot vs. time 400 
illustrating Thevenin voltage parameter definitions for an 
inverting CMOS gate. Shown are a gate input signal 4 02, a 
Thevenin voltage 404, a switching threshold delay (tO) 



406, and a slope delay [deltat] 408. The delay between 
the time when the gate input signal 402 reaches the logic 
switching threshold Vst and the time when the Thevenin 
voltage 404 begins to fall defines the switching threshold 

5 delay (to) 4 06. The delay between the time when the 

Thevenin voltage 404 begins to fall and the time when the 
Thevenin voltage 404 reaches zero defines the slope delay 
(deltat) 408. Similar Thevenin voltage parameters may be 
defined for other CMOS gates and input transitions. 

0 Using the Thevenin model, the output voltage Vo' 

of the CMOS gate 15 0 may be expressed as a function of 
time t by 

Vo'it) ^ y{t,tO,deltaf) (1) 

output transition 

t - to < 0 (2) 
0 < t - to < deltat 

t - to > deltat 

transition 

t - to < 0 (3) 

0 < t - to < deltat 
t - to > deltat 



Specifically, for a falling 

15 y{t, to, deltat) = Vdd 

y{t,tO,deltat) = Vdd - yO{t,tO, deltat) 

y{t,tO, deltat) = Vdd - yO{t,tO, deltat) 
+ yO{t - deltat, to, deltat) 



and for a rxsing output 
2 0 y{t,tO,deltat) = 0 

y{t,tO,deltat) = yO{t,tO, deltat) 

y{t,tO,deltat) = yO{t,tO, deltat) 
- yO{t - deltat,tO, deltat) 



where yO(u,v^,w) is the single positive ramp response for a 
25 capacitive load C and a driver resistance Rd, which may be 
found from the equation 



yO{u,v,w) = Vdd ■ 



Equations (1), (2), (3), and (4) are derived and 
explained in detail in the article "Performance 
Computation for Precharacterized CMOS Gates with RC 
Loads", Florentin Dartu, Noel Menezes, and Lawrence T. 
Pileggi, IEEE Transactions on Computer-aided Design of 
Integrated Circuits and Systems , Vol. 15, No. 5, May 1996. 

FIG. 5 is a voltage plot 500 illustrating 
typical output voltage response Vo^ftJ of the CMOS gate 
150 of FIG. 1 as a function of time. As the effective 
capacitance of the load 152 builds up charge, the output 
voltage Vo^ (t) increases from zero to 3 0 percent of the 
final voltage for a rising transition in a time t3 0 from 
the input switching threshold and continues increasing to 
50 percent of the final voltage in a time t50 measured 
from the input switching threshold. For a falling 
transition, the output voltage Vo' drops to 70 percent of 
the initial voltage in a time t70 from the input switching 
threshold. 

In contrast to methods that solve for the 
Thevenin parameters using a single effective capacitance, 
the following method uses two values of effective 
capacitance Ceffl and Ceff2. Ceffl is used to find either 
tSO (Ceffl) (for a rising transition) or t70 (Ceffl) (for a 
falling transition), and Ceff2 is used to find t50(Ceff2) 
to calculate the Thevenin parameters from the ramp 
response equations 



y(t50(Ceff2), tO, deltat) =0.5- Vdd 
y(t30(Ceff1), tO, deltat) =0.3 - Vdd (or 
y(t70(Ceff1), to, deltat) = 0.7 -Vdd ) 



(5) 



Initial estimates are made for Ceffl, Ceff2, tO, 
and deltat, for example: 



Ceffl = Ceff2 = .0001 pf (6) 
5 deltat = 5 ■ t50(Ceff2) - t30(Ceff1 ) 

to = t50(Ceff2) -0.69 Rd ■ Ceff2 - 0.5 ■ deltat 

Newton' s method may be used to solve the non- 
linear equations (5) from equations (1) - (4) for tO and 

10 deltat by using the estimates of Ceffl and Ceff2 and a 
delay lookup table according to well known techniques. 
The delay lookup table contains a range of delays 
calculated as a function of input ramptime and input 
capacitance and is contained in a precharacterized cell 

15 library created according to well known techniques. The 
calculated solutions for tO and deltat are substituted for 
the previous estimates of tO and deltat in the ramp 
response equations (5) until the solutions for tO and 
deltat converge to a desired accuracy. 

2 0 After computing the Thevenin parameters tO and 

deltat from the initial estimates, either the Foster model 
200 of FIG. 2 or the pi model 3 00 of FIG. 3 may be used 
according to well known techniques such as AWE (asymptotic 
waveform evaluation) to calculate a first delay, delayl, 

25 to t30 (Ceffl) or t70 (Ceffl) and a second delay, delay2 , to 
tSO (Ceff2) . 

The calculated values for delayl and delay2 are 
compared to corresponding delays delayl ' and delay2 ' found 
in the lookup table using the input ramptime and Ceffl and 
30 Ceff2, respectively. If delayl and delay2 have converged 
to within a desired accuracy of delayl' and delay2 ' , then 
the calculation of the Thevenin parameters tO and deltat 
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is complete. If not, then new values for Ceffl and Ceff2 
are found by a reverse interpolation of delayl and delay2 
from the delay lookup table, i.e., 

delay = f(input ramptime, capacitance) => ( 7 ) 

capacitance = g(input ramptime, delay) 

After obtaining new values for Ceffl and Ceff2, 
the Thevenin parameters tO and deltat are recomputed from 
the ramp response equations (5) , and new values for delayl 
and delay2 are calculated from the Foster or the pi model 
until delayl and delay2 converge to delayl' and delay2 ' 
within a desired accuracy, for example, 1 picosecond. 

FIG. 6 is a flowchart 600 for calculating the 
Thevenin parameters of switching threshold delay tO and 
slope delay deltat and effective load capacitances Ceffl 
and Ceff2 according to the method explained above . 

Step 602 is the entry point for the flowchart 

600 . 

In step 604, estimates for Ceffl, Ceff2, tO, and 
deltat are initialized. 

In step 606, the solutions to the following ramp 
response equations are calculated: 

y(t50(Ceff2), tO, deltat) =0.5- Vdd and either 
y(t30(Ceff1), tO, deltat) = 0.3 ■ Vdd or 

y(t70(Ceff1), to. deltat) = 0.7 -Vdd for tO and deltat. 

In step 608, the estimates of tO and deltat are 
compared with the calculated solutions. If the solutions 
for to and deltat have not yet converged to the estimates, 
then the estimates of to and deltat are replaced with the 
calculated solutions in step 610, and step 606 is 
repeated. 



If the solutions for tO and deltat have 
converged to the estimates, then the delayl is calculated 
as a function of t30(Ceffl) or t70(Ceffl) and delay2 is 
calculated as a function of t50(Ceff2) from a Foster or a 
pi model in step 612. 

In step 614, corresponding delays delayl' and 
delay2i for Ceffl and Ceff2 are found in the delay lookup 
table. 

In step 616, the calculated values for delayl 
and delay2 are compared to the corresponding delays 
delayl' and delay2 ' . If delayl and delay2 have converged 
to within a desired accuracy of delayl' and delay2 ' , then 
control transfers to step 622. 

If delayl and delay2 have not converged to 
within a desired accuracy of delayl' and delay2 ' , then new 
values for Ceffl and Ceff2 are found from the delay lookup 
table in step 618 . 

In step 62 0, the estimates of Ceffl and Ceff2 
are replaced with the new values and step 606 is repeated. 

Step 622 is the exit point for the flow chart 

600 . 

The flowchart 60 0 may also be embodied in a 
computer program product that includes a medium for 
embodying a computer program for input to a computer and a 
computer program embodied in the medium for causing the 
computer to perform at least one of the following 
functions : 

(a) initializing estimates of effective 
capacitances Ceffl and Ceff2, of a switching threshold 
delay to, and of a slope delay deltat; 



(b) solving ramp response equations for tO and 
deltat as a function of Ceffl and Ceff2; 

(c) comparing the estimates of tO and deltat 
with solutions for tO and deltat found in step (b) ; 

(d) replacing the estimates of tO and deltat 
with the solutions for tO and deltat if the solutions for 
to and deltat have not converged to the estimates of tO 
and deltat; 

(e) repeating steps (b) , (c) , and (d) until the 
solutions for tO and deltat converge to the estimates of 
to and deltat; 

(f) calculating a delayl as a function of 
t30 (Ceffl) or t70 (Ceffl) and a delay2 as a function of 
t50 (Ceff2) from a Foster or a pi model; 

(g) comparing delayl and delay2 to delays 
delayl ' and delay2 ' corresponding to Ceffl and Ceff2 in a 
delay lookup table; 

(h) finding new values for Ceffl and Ceff2 from 
a reverse lookup of delayl and delay2 in the delay lookup 
table if the calculated values for delayl and delay2 have 
not converged to delays delayl ' and delay2 ' ; 

(i) replacing the estimates of Ceffl and Ceff2 
in step (b) with the new values for Ceffl and Ceff2; and 

(j) repeating steps (b) through (i) until the 
calculated values for delayl and delay2 converge to delays 
delayl ' and delay2 ' . 

Other modifications, variations, and 
arrangements of the present invention may be made in 
accordance with the above teachings other than as 
specifically described to practice the invention within 
the spirit and scope defined by the following claims. 



